Universal Behavior in Heavy Electron Materials 
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We present our finding that an especially simple scaling expression describes the formation of a new 
state of quantum matter, the Kondo Fermi liquid (KL) in heavy electron materials. Emerging at T* 
as a result of the collective coherent hybridization of localized f electrons and conduction electrons, 
the KL possesses a non-Landau density of states varying as (1 — r/r*)^''^[l + ln(r*/r)]. We show 
that four independent experimental probes verify this scaling behavior and that for Celrlns the KL 
state density is in excellent agreement with the recent microscopic calculations of hybridization in 
this material by Shim, Haule, and Kotliar. 
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Scaling behavior in heavy electron materials was first 
suggested by Nakatsuji et al. [H, hereafter NPF, who 
found, through a careful analysis of experiments [2| on 
the thermal and magnetic behavior of the 115 materi- 
als, Cei_a;Laa;CoIn5 and its Ir counterpart, that these 
results could be interpreted using a two-fluid model of 
their specific heat and bulk susceptibility. One fluid, 
the emergent itinerant heavy electron component, was 
characterized by an order parameter, f(T/T*), that in- 
creased with decreasing temperature and scaled with the 
temperature, T*, at which itinerant behavior emerged; 
the second component was a spin liquid, characterized 
by an order parameter, [1 — /(T/T*)], and made up of 
weakly interacting local moments, that at low tempera- 
tures could be described by a collection of non-interacting 
Kondo impurities. The itinerant component possessed 
two striking properties: its spin susceptibility and specific 
heat were related by a temperature independent Wilson 
ratio, Rw, and both exhibited scaling behavior, varying 
logarithmically as T* /T. The NPF proposal was subse- 
quently shown to extend to a broad spectrum of heavy 
electron materials by Curro et al. hereafter CYSP, 
who showed that the NPF two-fluid model provided a 
natural explanation for the hyperfine anomaly found in 
materials for which NMR and ^SR measurements of the 
Knight shift (KS) do not track the measured bulk spin 
susceptibility. By interpreting its appearance as the onset 
of two-fluid behavior, their work provided an independent 
determination of T* and of the heavy electron spin sus- 
ceptibility that emerges below T* , although for CeCoIns, 
the NMR expression for this quantity differed somewhat 
from its NPF thermal determination. Taken together, 
the work of these authors suggests that in a Kondo lat- 
tice the hybridization of the localized f-electron magnetic 
moments with the conduction electrons becomes a collec- 
tive or global process at T* that leads to the formation 
of a new state of quantum matter, the Kondo Fermi liq- 
uid, or KL, whose density of states displays distinctly 
non-Landau scaling behavior between T* and a typically 
much lower cut-off temperature, Tg, below which another 



form of quantum order appears. As we shall see, it co- 
exists there with both the spin liquid and a "light" elec- 
tron Landau Fermi liquid made of those conduction elec- 
trons that do not hybridize with the local moments. 

In this letter, we present the results of an analysis that 
unites the thermal and magnetic determinations of the 
parameters that characterize KL behavior and leads to 
an especially simple scaling description of its density of 
states. We show that this KL density of states is seen in 
a number of other experimental probes of heavy electron 
behavior, and that quite remarkably, when our scaling ex- 
pression for the KL state density is applied to the heavy 
electron material Celrlus, it yields a state density that is 
in excellent agreement with the recent microscopic calcu- 
lations for this material by Shim et al. Q ■ The agreement 
provides a double validation — of the methods used in 
their microscopic approach and the results obtained here 
using a phenomenological approach — and we discuss its 
implications for heavy electron behavior above Tq, and 
the new quantum ordered states that emerge below Tq. 

Our approach to reconciling small differences in the 
details of the scaling results of NPF and CYSP for the 
Co 115 material, and obtaining a simple form for its KL 
density of states, is perhaps best explained by comparing 
the NPF two fluid description of the bulk susceptibility: 



X^f{T)xKL + [l-f{T)]xsL, 



(1) 



in which a Fermi liquid contribution from those "light" 
conduction electrons that do not hybridize is neglected, 
and the comparable expression for the Knight shift of a 
probe nucleus (apart from a constant offset -^"0)1 



K-Ko = Af{T)xKL + B[l - f{T)]xsL, 



(2) 



in which the spin liquid component is coupled to the 
probe nucleus by a transferred hyperfine coupling con- 
stant B, while the itinerant Kondo liquid component is 
coupled by a direct hyperfine coupling A. The emergent 
anomalous component of the Knight shift considered in 
CYSP can then be written as 



= K-Ko-Bx={A- B)f{T)xKL. 



(3) 
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FIG. 1: (Color online) Scaling behavior of the anomalous 
Knight shift for Ce-115 and some other heavy electron mate- 
rials (data from Refs. 0,01 and references in Ref. Q). The 
solid lines represent the KL density of states. 



On combining Eqs. ([T|) and ([3]), one gets a very simple 
formula for the order parameter f{T), 



f{T) = 1 - 



Xsl{T) 



(4) 



in which a = \/{A~ B) and the only unknown quantity 
is XSL- The behavior of xsl is known in two limits: at 
high temperatures, T > T* , xsl is given by the Curie- 
Weiss law, xcwiT) ^ {T + r*)~^; at low temperatures, 
T < T*/4, say, because the local f-moments are strongly 
screened by the weakly hybridized conduction electrons, 
it reduces to a collection of non-interacting single Kondo 
impurities, an assumption born out by the NPF results. 
A simple linear interpolation between these limits then 
yields xsl, while a is determined by comparing the re- 
sults obtained with above approach with those found by 
a self-consistent determination of f(T) based on the NPF 
approach. The result we obtain for the order parameter. 



/(T/T*) = /(O) ( 1 - ^ 
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(5) 



can be combined with our result for the KL quasiparti- 
cle effective mass determined from the specific heat and 
Knight shift fits. 



'■KL 



= m,, ( 1 + In — 



(6) 



(1-T/T*) [1+ln(T*/T)] 




(a) 



(1-T/T*) [1+ln(T*/T)] 



(b) 



Ce,.,La,Coln5 



1.2 



FIG. 2: (Color online) Scaling behavior of the anomalous Hall 
coefficient for (a) Ce-115 materials and (b) Cei-2,Laa;CoIn5. 
The solid lines are the KL density of states. 



seen in xk l , where mu is the quasiparticle effective mass 
at T* , to yield the reduced density of states for the Kondo 
liquid. 
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Fig. [T] shows that our result provides an excellent fit to 
the values of Kanom obtained from NMR experiments for 
both the 115 family of heavy electron materials and all 
others that display a Knight shift anomaly. 

Three recent experiments on the 115 materials provide 
additional evidence for the presence of the Kondo liquid 
and a description of its density of states using Eq. ([71): the 
anomalous Hall effect, tunneling experiment, and Raman 
scattering. 

Anomalous Hall effect 

Hundley et al. [7| have found that the Hall effect in 
the Ce-115 compounds does not follow the standard skew 
scattering theory, based on single Kondo impurity scat- 
tering, that explains the behavior of many heavy electron 
materials [8|. In this theory, the influence of applied ex- 
ternal fields on the local Kondo impurities produces a 
left-right asymmetry that is seen as a skew scattering 
contribution to the Hall effect. To understand the de- 
parture of their results from this theory, Hundley et al. 
introduce the dimcnsionless function. 



a{T) = (i?g'=(T) - i?f '^"'(T)) /i?]^°(T) 



(8) 
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FIG. 3: (Color online) Scaling behavior of the conductance 
asymmetry in the tunneling experiment on CeCoIns . The KS 
anomaly is also shown for comparison. The solid line indicates 
the KL density of states. 



to describe additional effects originated from the local- 
ized Ce ions. Here Rjj'^{T) is the Hall coefficient of the 
Ce-115 compound and Rjf{T) is the Hall coefficient of 
the corresponding La compound. Hundley et al. find that 
the skew scattering term Rj^^'^{T) makes a minor con- 
tribution to the total Hall effect and that while a{T) is 
nearly constant at high temperatures, it increases sharply 
below a characteristic onset temperature that we identify 
with T*. On analysing their results, we find that a{T) 
exhibits exactly the following scaling behavior: 



v{T)-a{300K)^ Pkl{T/T*). 



(9) 



As may be seen in Fig. [D^a), the scaled a{T) — a{300K) 
for all three Ce-115 materials falls upon the universal 
density of states of the Kondo liquid when one takes 
the onset temperature T* for the anomalous Hall effect 
to be 53 K for CeColns, 31 K for Celrlng, and 18 K for 
CeRhlns, results that are comparable to the values of 
T* for these materials obtained from an analysis of their 
magnetic and thermal behavior. As shown in Fig. [^I^b), 
similar scaling behavior is also found in Cei_2,Laa;CoIn5 
0. In these materials, the onset temperature of the Hall 
anomaly again agrees very well with the earlier results 
from the magnetic susceptibility and the specific heat and 
displays clearly the linear decrease of T* with increasing 
La doping that was found by NPF. We conclude that 
Hall measurements not only provide a direct measure of 
the KL state density, but offer a reliable independent es- 
timate of T*. 

Tunneling experiment 

Kondo liquid behavior has also been found to play a 
role in tunneling conductance. In the recent experiment 
by Park et al. on CeCoIns [lo| . a conductance asymme- 
try, defined as the ratio of the differential conductance 
at -2 mV and 2 mV, is found to develop at the charac- 
teristic KL temperature T* and to increase with decreas- 
ing temperature down to the superconducting transition 
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FIG. 4: (Color online) Fano asymmetry parameter \q\~^ from 
the Raman spectra of CeCoIns compared to the KL density 
of states. 



temperature Tc (2.3 K). As Curro [ll| has pointed out, 
their measured asymmetry follows the KS anomaly mea- 
sured by CYSP and as may be seen in Fig. [31 both follow 
the universal KL density of states for temperatures above 
the cut-off temperature. To ~ 10 K, that marks the end 
of Kondo scaling behavior, while their mutual agreement 
continues down to even lower temperature. This tells us 
that the breakdown of universality does not necessarily 
indicate the breakdown of the two fluid description, but 
signals instead a transition of the KL state to some other 
low temperature quantum state. 
Raman scattering 

Very recently, Martinho et al. [l3| have used Raman 
scattering to probe the behavior of CeCoIn5. They note 
that the Fano line shape of the Raman spectral function 
probes the coupling of phonons to itinerant electrons and 
that the Fano asymmetry parameter q becomes markedly 
temperature dependent below a temperature of order T* . 
This led us to compare it with our KL state density, with 
the results shown in Fig. [H While not as convincing 
as the anomalous Hall effect or tunneling conductance 
asymmetry, these results suggest that in the future Ra- 
man scattering might offer a useful probe of Kondo liquid 
behavior. 

At this stage it is natural to inquire what, if any, theo- 
retical basis exists for our phenomenological description 
of Kondo liquid behavior. The answer comes in the very 
recent work of Shim et al. hereafter SHK, who carry 
out a first principles calculation of the localized to itin- 
erant transition in Celrlns using Dynamical Mean Field 
Theory in combination with Local Density Approxima- 
tion (LDA-I-DMFT). In Fig. [5] we compare their calcu- 
lated density of states at the quasiparticle peak with our 
universal KL density of states; T* is chosen to be 31 K, 
as suggested by the anomalous Hall data. Since in the 
SHK result there is no clear onset temperature at which 
the density of states starts to increase, we simply com- 
pare the low temperature (T < T*) portions and find a 
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FIG. 5: (Color online) LDA+DMFT density of states of 
the quasiparticles in Celrlns compared to the KL density of 
states. We take the onset temperature T* = 31 K from the 
Hall anomaly. 



linear relation between them. 

The excellent agreement between the two results indi- 
cates that SHK have captured the Kondo liquid state 
density without having to resort to a cluster calcula- 
tion below T* , a truly remarkable finding that provides 
quite strong evidence for the validity of both their mi- 
croscopic approach to hybridization and our phenomeno- 
logical scaling expression for the KL state density. Since 
there is no experimental evidence for a change in the spe- 
cific heat at the hyrbridization phase transition (i.e. noth- 
ing special happens at T*), and since local hybridization 
can influence the density of states up to very high tem- 
peratures, it is to be expected that an exact calculation 
of the state density would do just what their calculation 
does, i.e. go smoothly through T* and yield our scaling 
results below T*. 

It is useful to make a distinction between the local 
(i.e. single Kondo impurity) hybridization that must be- 
gin at temperatures well above T* and the universal 
global hybridization behavior that begins at T*. The 
local hybridization is expected to be highly anisotropic 
and complex because the Kondo coupling, J, of the sin- 
gle Kondo impurities to different parts of the Fermi sur- 
face will be different; this is just what SHK flnd. T* 
is the temperature at which, as a result of feedback ef- 
fects on the coupling between individual hybridization 
channels, hybridization becomes a "simple" collective co- 
herent universal phenomenon, characterized by an or- 
der parameter, /(T/T*), that goes as (1 - T/T*)3/2, 
and a logarithmic increase in the average effective mass, 
to|^^ = m;i[l -|- ln(T*/T)], where rrih incorporates the 
hybridization occurring above T*. We further note that 
NPF argue that the physical origin of the energy scale 
T* that characterizes KL behavior is the nearest neigh- 
bor coupling between the localized electrons, a view that 
is supported by its appearance as the Curie- Weiss param- 



eter in the high temperature susceptibility for CeCoIns 
and other materials in which crystal field effects play out 
at substantially higher temperatures. 

We conclude that once one has identified T* for a given 
heavy electron material, the interacting components nec- 
essary to understand its transport behavior at all temper- 
atures as well as the appearance of quantum order below 
Tq can now be specified. There will in general be three: 
a heavy electron non-Landau Kondo Fermi liquid whose 
quasiparticles have the temperature dependent effective 
mass given by Eq. ([6]); a "light" electron Landau Fermi 
liquid corresponding to electrons on those parts of the 
expanded Fermi surface that have not hybridized individ- 
ually or collectively; and a residual spin liquid made up 
of the screened interacting local moments. So to under- 
stand the behavior of heavy electron materials between 
T* and Tq one needs first of all to combine our scaling de- 
scription of the Kondo liquid behavior and the two fluid 
model to follow the temperature evolution of all three 
components, while it is their mutual interaction, in the 
presence of global hybridization, that determines trans- 
port properties above Tq, and the nature of the ordered 
state below Tq. 
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